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Abstract. We discuss the generic chiaracteristics of stocliastic particle acceleration by a fully developed turbulence spectrum 
and show that resonant interactions of particles with high speed waves dominate the acceleration process. To produce the 
relativistic electrons inferred from the broadband spectrum of a few well-observed shell-type supernova remnants in the 
leptonic scenario for the TeV emission, fast mode waves must be excited effectively in the downstream and dominate the 
turbulence in the subsonic phase. Strong collisionless non-relativistic astrophysical shocks are studied with the assumption 
of a constant Aflven speed. The energy density of non-thermal electrons is found to be comparable to that of the magnetic 
field. With reasonable parameters, the model explains observations of shell-type supernova remnants. More detailed studies 
are warranted to better understand the nature of supernova shocks. 
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TURBULENCE CASCADE AND 
STOCHASTIC PARTICLE 
ACCELERATION 

In the Kolmogorov phenomenology, the free energy 
dissipation rate is given hy Q = Cipu^/L, where Ci 
is a dimensionless constant, and u and L are the eddy 
speed and the turbulence generation scale, respectively. 
The eddy turnover speed and time at smaller scales are 
given respectively by v^j^j{k) = 47zW{k)k^ <^ k^^^^ 
and x,d,i{k) = 2n/kvedd = n'l^{Wk^)-'l^ - 
k-^l^ , where W{k) = {u^ /An){2n/Lfl^k-^^l^ = 
{'Xn)-'^{2nQ/C\pfl^k-^^l^ ^ /t^''/^ is the isotropic 
turbulence power spectrum, k — 2n/l is the wave num- 
ber and I is the eddy size. From the three-dimensional 
Kolmogorov constant C ~ 1.62 ifioll . we obtain 
Ci — 2;r/C^/^ = 3.05. At the turbulence generation 
scale k,„ = 27t/L, v,,ui ^ u, Q ^ 2Cip[47:W^k^^]^/^ = 
^iP'^'lddi^)/'^edd{k) , and the total turbulence energy 
is given by JW{k)47i:k^dk = {3/2)u^. The turbulence 
decay time is therefore given by T^/ = 3t„/^(A;,„)/Ci , i.e., 
eddies decay after making 3Cf ' ^ 1 turns. 

We are interested in the acceleration of particles 
through scattering randomly with heavy scattering cen- 
ters with the corresponding acceleration time t^c — 
V[3v^/v^^^(^)] , where Tjy- — 2k /kv ~ l/v is the scat- 
tering time, V is the particle speed, and we have as- 
sumed that the scattering mean free-path is equal to I. 
For the above isotropic Kolmogorov turbulence spec- 
trum, Tac{k) = 3v/2Wk^ fc^'/^. To have significant 
stochastic particle acceleration (SA), the acceleration 
time Tac{k) should be shorter than the turbulence de- 



cay time T^/, which implies u^ > C\vvgiiii{k) . So, in gen- 
eral, the SA is more efficient at smaller scales. The on- 
set scale of the SA is given by kc = (Ci v/m)^^,,, . There- 
fore, to produce energetic particles with a speed of v by 
a Kolmogorov spectrum of scattering centers, the turbu- 
lence must have a dynamical range greater than Dj^ = 
{Civ/u)\ 

In the Kraichnan phenomenology, the turbulence de- 
cay is suppressed by the wave propagating effect with 

Q = CipuyLvF, W{k) = {u^l\n)k]Jh-^l^, Tedd - 
k^^^^ ,Vi,dd °^ k^^l"' ,Tac °^ A;^'/^, and the turbulence de- 
cay time Td = 3Tedd{km)vF /Ciu , where the wave speed 
Vf ^ u. To have significant acceleration through scat- 
tering with the eddies, the dynamical range of the tur- 
bulence must be greater than D^: = (Civ/vf)^ , which 
is much less than = {Civ/u)^. The resonant inter- 
actions of particles with waves can be much more effi- 
cient in accelerating particles in this case. For a wave 
speed Vf independent of k, the acceleration time is given 
by Tac = 3TscV^/vjr °= . To have significant accelera- 
tion, the scattering mean free path of the particles must 
be shorter than (v^ /Ci vu'^)L = oj^^D^^^^^L . 

Several Shell-Type Supernova Remnants (STSNRs) 
have been observed extensively in the radio. X-ray, and 
TeV bands. X-ray observations with Chandra, XMM- 
Newton, and Suzaku and TeV observations with HESS 
have made several surprising discoveries that challenge 
the classical diffusive shock particle acceleration model 
iS H. The SNR RX J1713.7-3946 is about t 1600 
years old ^ with a radius of ~ 10 pc and a distance 
of D ~ 1 kpc. By fitting its broadband spectrum with 
an electron distribution of / °^ 7^''exp — (7/7^)'/^, we 



find that p = 1 .85, B = 12.0 ^G, /. m^c^ = 3.68 TeV, and 
the total energy of relativistic electrons with the Lorentz 
factor 7 > 1800 Ee = 3.92 x 10"^^ erg (Fig.©. 

The X-ray emitting electrons have a gyro-radius of 
rg ~ 10^^ cm, which shouldn't be shorter than the scatter- 
ing mean free path. To produce these electrons through 
the SA, the turbulence must be generated on scales 

3/2 —2/3 

greater than Di^rg, D^rg, and D„ D^. Vg for the non- 
resonance Kolmogorov, Kraichnan phenomenology, and 
the resonant interactions, respectively. For STSNRs, u ^ 
Vf ^ 0.01c, Dj^rg ^ 10 kpc, which is much larger than 
the radii of the remnants. The SA by eddies with a Kol- 
mogorov spectrum is therefore insignificant. D^rg ^ 30 

3/2 —2/3 

pc, which is also too thick. Z)vJ Z)^ ' ^ 0.1 pc, which 
is much greater than the particle inertial length and may 
be generated through the Kelvin-Helmholtz instabili- 
ties or cosmic ray drifting upstream (i3,|5(]. Therefore if 
relativistic electrons from the STSNRs are accelerated 
through the S A, they must be energized through resonant 
interactions with high speed plasma waves. Low speed 
waves also require a large turbulence dynamical range to 
accelerate particles. 



SHOCK STRUCTURE, WAVE DAMPING, 
AND STOCHASTIC ELECTRON 
ACCELERATION BY FAST MODE 
WAVES IN THE DOWNSTREAM 

We next study the SA in the shock downstream by 
weakly magnetized turbulence with the Alfven speed 
va = {B^ /A-TtpY^^ ^ where B, and p are the magnetic 
field, and mass density, respectively. For strong non- 
relativistic shocks with the shock frame upstream speed 
U much higher than the speed of the parallel propagating 
fast mode waves in the upstream Vf — (v^ + 5v|/3)'/^, 
where v| = P/p is the isothermal sound speed and P is 
the gas pressure, mass, momentum, and energy conser- 
vation across the shock front require 



(1) 



where we have assumed that the turbulence behaves as an 
ideal gas and ignored the wave propagation effects. The 
shock structure can be complicated due to the present of 
turbulence. We assume that the turbulence is isotropic 
and has a generation scale of L, which does not change 
in the downstream. The speeds vj, v^, and u therefore 
should be considered as averaged quantities on the scale 
L. va depends on the upstream conditions and/or the 
dynamo process of magnetic field amplification|[ll IH. 
Here we assume it a constant in the downstream. One 
can then quantify the evolution of other speeds in the 
downstream. 
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FIGURE 1. Evolution of the eddy speed u and speed of 
parallel propagating fast mode waves vp in the downstream for 
VA = 0.0633^/. 



For the Kolmogorov phenomenology. 
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Near the shock front, we denote the isothermal sound 
speed and Aflven speed by vsq and v^o, respectively. 
Then the eddy speed at the shock front is given by 
fli/2{//4 with fl = 3 - 16v|o/t/2 - 
equation (|2| from the shock front i 
(x > 0), we then have 



32vIq/5U^. Integrate 
X = 0) to downstream 
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As mentioned in the previous section, to produce the ob- 
served X-ray emitting electrons in the STSNRs through 
the SA processes, fast mode waves must be excited ef- 
ficiently. The MHD wave period is given by Tf (fc) ~ 
In/vpk. Then the transition scale from the Kolmogorov 
to Kraichnan phenomenology occurs at Tf {k, ) = t„m (kt ) 
or VF — Vg^^(fef)[2J. We then have 



kt = {u/vpfkn 



(6) 



For k > k, > k,„, the turbulence spectrum in the iner- 
tial range is given by W{k) = u^{4n)^^km^k, '^^^^^/^ — 
(4;r)"^vy^M-'/2^^/2^-7/2 ^ Although the turbulence en- 
ergy exceeds (3/2)m2 when the wave propagation effect 



is considered, we still assume that the enthalpy of the tur- 
bulence is given by (5/2)m2 for < M SO that equation 
([T]i and the above solutions for the speed profiles remain 
valid. 

In the subsonic phase with Vf > u, we assume that fast 
mode waves can still be excited efficiently to maximize 
the SA efficiency. Then the Kraichnan phenomenology 
prevails and 



w{k) = u^i47t)-hii^k-y^, 

3Udu{xf _ Cim(jc)'* 
8d^ ^ 
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where from equation ([TJ one has vp = 

[5U^/16 + vl/3-5u^{x)/3Y^^ . These equations 
can be solved numerically to get the speed profiles 
in the subsonic phase. Figure [T] shows the vp and u 
profiles with va = v^o = 0.0633t/ in the downstream 
and V50 = VAO < U. 
In summary, 
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in the Kraichnan regime. The coUisionless damping 
starts at the coherent length of the magnetic field l^/ = 
In/k^, where the period of Alfven waves In/kvA 
is comparable to the eddy turnover time, i.e., v? = 



AnW{kd)k] 
have 



mm we 

k^ = [u^mm{vf,u)/vA]k,„. (10) 

For a fully ionized hydrogen plasma with isotropic par- 
ticle distributions, which is reasonable in the absence of 
strong large scale magnetic fields, the transit-time damp- 
ing (TTD) rate is given by 10, 01 
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where k^, Tg, Tp, nig, nip, 9, CO, and = kcos9 are 
the Boltzmann constant, electron and proton tempera- 
tures, masses, angle between the wave propagation di- 
rection and mean magnetic field, wave frequency, and 
parallel component of the wave vector, respectively. 
The first and second terms in the brackets on the right 
hand side correspond to damping by electrons and pro- 
tons, respectively. For weakly magnetized plasma with 
Va < V5, proton heating always dominates the TTD for 



k^Tp/mp. If Va does not change dra- 



matically in the downstream, the continuous heating of 
background particles through the TTD processes makes 
Tp {mp/mg)Tg since the heating rates are proportional 




FIGURE 2. The turbulence spectra v^^j{k) at several loca- 
tions in the downstream as indicated. The Kolmogorov, Kraich- 
nan, and damping ranges are indicated for the supersonic phase 
spectrum with x = O.OIOL. At the other locations, the turbu- 
lence is subsonic and there are only Kraichnan and damping 
ranges. The sharp drops of the turbulence spectra in the damp- 
ing range are due to the onset of thermal damping at the coher- 
ent length of the magnetic field iTt/kii. 



to (mr)'/^, where m and T represent the mass and tem- 
perature of the particles, respectively. We see that par- 
allel propagating waves (with sin0 ~ 0) are not sub- 
ject to the TTD processes and can accelerate some par- 
ticles to very high energy through cyclotron resonances. 
Obliquely propagating waves are damped efficiently by 
the background particles. Although the damping rates for 
waves propagating nearly perpendicular to the magnetic 
field (cos ~ 0) are also low, these waves are subject 
to damping by magnetic field wandering |6]. The turbu- 
lence power spectrum cuts off sharply when the damping 
rate becomes comparable to the turbulence cascade rate 



r = ^eddli^F + ^edd) - %dd^F 12]. One can define a 
critical propagation angle 9c{k), where Aj- ( 0^ , ^) = r(^)- 
Equations (|9]l and ( fTTT l give 



^Ai^d 
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l^lH'^l^sVFk^l^ - COS0, 



exp 
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where the electron damping term has been ignored. The 
turbulence spectra at several locations in the downstream 
are shown in Figure |2] 

The escape time of relativistic particles with v ~ c, 
where c is the speed of light, from the particle acceler- 
ation region is given by Xg^c = (i-^/4c^)T(c. and the spec- 
tral index of the accelerated particles in the steady state 
is given by 
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FIGURE 3. Evolution of the acceleration efficiency rj (dot- 
ted), cutoff Lorentz factor Yc (dotted-dashed), spectral index p 
(dashed), and T = tac/tj (solid) in the downstream for = 
0.0633f/ and U = 0.01c. The particle acceleration is signifi- 
cant for T < 1, i.e., between the two vertical dashed lines in- 
dicating xi and X2. For Yc, we have assumed that B = 12/xG 
and L = 0.15pc. For 7], the two thin lines are for V = 0.015c 
(higher) and 0.0067c (lower). See the following section for de- 
tails 
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We note that for independent of x in the downstream, 
p reaches its minimum at the transonic point xq, where 
Vf ~ u. The maximum energy that particles can reach 
though resonant interactions with these parallel propa- 
gating waves is given by 



InqB _ 
nieC^kj mgC^u^ mm(vF ,u) 
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where q is the elementary charge units. The ratio of the 
dissipated energy carried by non-thermal particles to that 
of the thermal particles should be greater than 
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where e — 2.72, since the isotropic turbulence 
with k < kj can also accelerate particles with the 
Lorentz factor 7 > jc- To have efficient accelera- 
tion of relativistic particles, the turbulence decay 
time Trf = 3max(M, )L/CiM^ should be longer 
than the acceleration time Xac — (3c^/v^)Tic = 
6nc/vpkd = 3cviL/vpU^min{vF,u) , which implies 
max(M, Vf )L/CiM > cv\L/vjfU^mm{vf,u), i.e., 
Ci < M^v^ /cv\ . There are at most two locations 
xi < X2 in the down stream, where T = Tac/Tj = 1 
and Ci = M^v^/cv^. In combination with equa- 
tion ( fTSl ), significant particle acceleration occurs for 



FIGURE 4. Normalized nonthermal electron distribution / 
produced at several locations in the downstream. 




FIGURE 5. The distributions of nonthermal electrons F in 
the downstream. 



< [9/4 + 1 2 max (m2 , ) /m^C^] ' - 1 /2 . The par- 
ticle acceleration in the supersonic phase, i.e., Vf < u, 
therefore produces very hard electron distributions with 
p < 1.39 for Ci — 3. Softer electron distributions have 
to be produced in the subsonic phase. Figure |3] shows 
the evolution of rj, = '^^ac/td in the down- 

stream for U = 0.01c. The profiles of vp/U and u/U 
only depend on va/U. So is the profile of 77. T and p 
also depend on the absolute value of U. To obtain 7, 
one needs to know L and B as well. In the extremely 
supersonic phase with Vf <C u, the SA is negligible. The 
SA is significant only after the plasma is already heated 
up so that VF ^ u. In the late subsonic phase, u ^ vf, the 
SA is also insignificant since most of the free energy of 
the system has been converted into heat. 



The particle distribution may be approximated rea- 
sonably well with /(x,7) o<: 7^^'^'exp-[7/7c(x)]'/^|[3t]. 
Then the distribution of non-thermal particle in the 
downstream 

F{x,Y)^ r f{x\Y)riix'){4Q/m,c^U)dx' (16) 

where X,"^/„,^ 7/(-^', 7)d7 = 1, and ymeC^F{x, 7)d7 

gives the energy density of non-thermal particles at x. If 
/cv\ < C\ at the sonic point .icq, thenjco < -^i and there 
will be no particle acceleration in the supersonic phase. 
Figures|4]and|5]show the normalized electron distribution 
/ and F at several locations in the downstream, respec- 
tively. 




RESULTS 

Here, we use the SNR RX J17 13.7-3946 as an example 
to demonstrate how the SA by fast mode waves accounts 
for the observed broadband spectrum. Figure |6] shows the 
best fit with VA/f/ = 0.0633, L = 4.71 x 10'^cm,B= 12 
jUG and U ~ 0.01c. Comparing to the thin dashed line, 
which is derived by assuming an electron distribution 
°^ Y^P exp — (7/ 7c) there is a radio spectral bump due 
to electron acceleration relatively far from the shock 
front (see Figs. |4] and |5]l. In our model, there are five 
parameters: B, U, va, L, and the equivalent volume of 
a uniform emission range. The last is 4 times bigger than 
the volume of the SNR, suggesting higher nonthermal 
electron densities in the interior of the remnant than near 
the shock front. The observed radio to X-ray spectral 
index. X-ray to TeV flux ratio, location of the X-ray 
cutoff, and bolometric luminosity of the source give four 
constraints, which leads to one more degree of freedom. 
Our model fit to the spectrum therefore is not unique. 
However, B is uniquely determined by the ratio of the 
X-ray to TeV flux. To reproduce the observed spectral 
shape, the profiles of p, 7., and rj should not change, 
which implies that v\c^ m"* and L °^ /v\ at the 
transonic point xq. For va <C f/, m is proportional to 
U . We therefore obtain va and L as functions of U 
as indicated in Figure IT] The density can be derived 
from B and v^, and the overall acceleration efficiency 
is defined as rf^q = .[^i iQ*^/ lo Q'^- Nearly identical 
spectrum can be obtained for parameters on these lines. 
We note rj oc v^/v^Vf °^ v^/U^ °^ U^^^. The acceleration 
is more efficient in the earlier phase of the remnant 
evolution. The two thin dotted lines in Figure[3]show the 
dependence of T] on f/. 
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